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ABSTRACT 

Recently, Tseliakhovich and Hirata (2010) showed that during the cosmic Dark Ages the baryons 
were typically moving supersonically with respect to the dark matter with a spatially variable Mach 
number. Such supersonic motion may source shocks that inhomogeneously heat the Universe. This 
motion may also suppress star formation in the first halos. Even a small amount of coupling of 
the 21cni signal to this motion has the potential to vastly enhance the 21cm brightness temperature 
fluctuations at 15 < z < 40 as well as to imprint distinctive acoustic oscillations in this signal. 
We present estimates for the size of this coupling, which we calibrate with a suite of cosmological 
simulations of the high-redshift Universe using the GADGET and Enzo codes. Our simulations, 
discussed in detail in a companion paper, are initialized to self-consistently account for gas pressure 
and the dark matter-baryon relative velocity, Vhc (in contrast to prior simulations). We find that the 
supersonic velocity difference dramatically suppresses structure formation on 10 — 100 comoving kpc 
scales, it sources shocks throughout the Universe, and it impacts the accretion of gas onto the first 
star- forming minihalos (even for halo masses as large as 10^ Mq). However, prior to reheating by 
astrophysical sources, we find that the Vbc~sourced temperature fluctuations can contribute only as 
much as « 10% of the fluctuations in the 21cm signal. We do flnd that Ubc in certain scenarios 
could source an 0(1) component in the power spectrum of the 21cm background on observable scales 
via the X-ray (but not ultraviolet) backgrounds produced once the first stars formed. In a scenario 
in which 10^ Mq minihalos reheated the Universe via their X-ray backgrounds, we find that the 
pre-reionization 21cm signal would be larger than previously anticipated and exhibit more significant 
acoustic features. Such features would be a direct probe of the first stars and black holes. In addition, 
we show that structure formation shocks are unable to heat the Universe sufficiently to erase a strong 
21cm absorption trough at 2; ^ 20 that is found in most models of the sky-averaged 21cm intensity. 
Subject headings: cosmology: theory — first stars — galaxies: high redshift - stars: Population HI - 
galaxies: formation 



1. INTRODUCTION 

Sometime between the redshifts of 30 and 10, it is 
thought that the Universe transitioned from a pristine 
landscape, where the distribution of matter was calcu- 
lable from the cosmological initial conditions alone, to 
a vastly more complex system in which stars abound. 
Once these stars formed, within just a few hundred 
million years they had reheated and reionized the en- 
tire Universe. These cosmic times can be observed di- 
rectly via redshifted 21cm radiation from neutral hydro- 
gen. In fact, predictions are that the 21cm line is vis- 
ible in absorption against the Cosmic Microwave Back- 
ground (CMB) from 15 < z < 40, with a larger ampli- 
tude for its mean brightness ter nperature than from any 
other cosmological ep och (e.g., iPritchard fc Loebl 120081 : 
iFurlanetto et al.l [20061 ). Motivated by these predictions, 
the radio instruments LEDA, DARE, and LOFAR are 
being designed to observe 21cm emission from the era 
rior to reiq i iization, the so-called cosrn i c Dark Ages 
Burns et al.l 120111: iBernardi et"all l2012t iHarker et all 
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l2OT0l ).^ These efforts are in addition to those aiming 
to detect the 21cm signal from the Reionization Epoch 
(targeting z ~ 6 - 12), which includes EDGES, PAPER, 
MWA, LOFAR (with t heir high band antennae ) , and 
GMRT (iBo wman & R ogeTi [2010t iParsons et al.l l2010l : 
IHarker etal.n201(k iPaciga et al.l l2011h. EDGES. LEDA. 
and DARE are targeting the sky-averaged 21cm signal, 
whereas the other instruments are attempting to measure 
fiuctuations in the 21cm intensity on ^ 10 comoving Mpc 
scales. Plans for the next generation of instruments are 
also in development.^ 

Yet, despite the prospects for observations of 21cm 
emission from 10 < z < 40, key theoretical questions 
regarding these times remain unanswered - even ques- 
tions which are independent of the large uncertainties 
inherent to modeling the astrophysical processes that 
impact the 21cm signal. For example, it is unclear 
whether weak structure formatio n shocks would have sig- 
nificantly heated the cosm ic gas (IGnedin fc Sha ver 200| 
IFurlanetto fc Loebl I2004D . In fact. iGnedhTfc Shaverl 
( 200l " found that the reheating from structure forma- 
tion shocks significantly impacted their simulated 21cm 
absorption signal. The sky-averaged 21cm brightness 
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temperature, which LEDA and DARE aim to observe, 
is hkely to be inversely proportional to the gas tempera- 
ture at the time of emission, so that colder temperatures 
would result in a larger signal. Here we show that struc- 
ture formation shocks are not likely to suppress this ab- 
sorption signal, in con trast to Gnedin fc Shaver 1200A) . 

In addition, recentlv lTseliakhovich fc Hiratal ( 201Clt ) ar- 
gued that a previously unconsidered effect could impress 
large spatial fluctuations in the 21cm signal: At the time 
of recombination, the cosmic gas was moving with re- 
spect to the dark matter at a root mean square velocity 
of 30 km s~^ and in a coherent manner on < 10 comov- 
ing Mpc separations. These initial velocity differences 
translate into the dark matter moving through the gas 
with a typical Mach number of « 2 from z = 150 
until the time the gas was reheated. Semi-analytic mod- 
els predict that astrophysical backgrounds reheat ed the 
Universe somewhere between z — 10 and 20 ( Furlanettol 
l200a IPritchard fc Loebll20Tol) . 

The baryon-dark matter supersonic motion 
could suppress the f or niation of the firs t stars 
(iTsehakhovich fc Hirat al 12011 IMaio et all I20TTI: 
l^acy et al. 2011; Na oz et al.l 1201 111 It could also 
introduce spatially variable kinetic temperatures by 
the entropy generated in shocks that these supersonic 
flows sourced. Both of these effects would result in 
fluctuations in the 21cm signal. In fact, the 21cm 
correlation function may have a significantly larger 
amplitude on the 10 — 100 comoving Mpc scales that 
observations are most sensitive to if this differential 
velocity impacts the 21cm background ([Dalai et al.l 
I2OIOI) . Already, many 21cm models tha t do not account 
for the ITsehakhovich fc Hiratal (|2010[ ) effect find that 
the fluctuations in the 21cm background from z = 20 
are likely to be no more difhcu lt to detect than those 
from z = 10 (|Pritchard fc Loebl l2008). Existing interfer- 
ometric efforts to detect highly redshifted 21cm emission 
are primarily focused on z = 6 — 10. If the signal from 
z = 20 were further boosted by the dark matter-baryon 
velocity differential (a scenario the results of this paper 
support under specific circumstances), then more effort 
should be channeled towards detecting 21cm fluctua- 
tions from this earlier epoch. Furthermore, the power 
spectrum of the differential velocity features significant 
acoustic features similar to those in the CMB. These 
features could provide a more distinctive signature that 
would aid the separation of the 21cm signal from the 
foregrounds. 

The aim of this paper is to estimate the coupling 
strength of the 21cm signal to the dark matter-baryon 
velocity differential. To do so, this study presents semi- 
analytic estimates for the size of this coupling, which 
are calibrated with numerical simulations of the Dark 
Ages. Our si mulations, described in de tail in a com- 
panion paper (jO'Leary fc McQuinnll2012l Paper I), em- 
ploy more physical cosmological initial conditions than 
the simulations used in prior studies. Unlike previous 
studies, our simulations are initialized with a transfer 
function that consistently incorporates the dark matter- 
baryon differential velocity. (On the scale of our sim- 
ulations, < 1 Mpc//i, this differential velocity is a uni- 
form wind to an excellent approximation.) In Paper I, 
we show that these improvements provide a much better 
match to the linear evolution of cosmological perturba- 



tions. In addition, we run both the Enzo (|0'Shea et al.l 
I2OOI and GADGET (|Springel et al.ll200l cosmological 
codes, with many different box sizes and particle num- 
bers/grid sizes, to explore the robustness of our results. 

This paper is organized as follows: Section |2] pro- 
vides a short introduction into the 21cm signal, its 
detectability, and introduces a parameterization for 
how this signal could couple to the dark matter-baryon 
velocity differential. We then provide estimates for 
how the 21cm signal is impacted by this velocity 
differential under various assumptions, using a suite 
of simulations to calibrate these assumptions (Section 
This study assumes a flat ACDM cosmological 
model with = 0.27, = 0.73, h = 0.71, erg = 0.8, 
Hs = 0.96, Yho = 0.24, and Qh = 046, consistent 
with recent measurements (jLarson "ern]l201lh . We wiU 
subsequently abbreviate proper Mpc as pMpc and Mpc 
and kpc will be reserved for comoving lengths. Some of 
our calculations use the Sheth-Tormen mass function, 
for which we adopt the paramete rs p — 0.3, a — 0.75, 
A 0.322 (ISheth fc Tormenll200l . 



Wh ile this study was nearing completion, IVisbal et al.l 
(|2012l ) p ublished a related study using semi-numeric 
methods. IVisbal et al.l (|2012| ) investigated one model for 
how the baryon-dark matter velocity difference could im- 
pact the 21cm signal. We provide here a more complete 
census of the different scenarios where the relative veloc- 
ity could contribute to the 21cm signal. 

2. THE PRE-REIONIZATION 21CM SIGNAL 

This section provides a brief introduction to the tem- 
poral evolution of the redshifted 21cm signal. This in- 
troduction differs slightly from prior e xpositions (e.g., 
iFurlanettol [20061 : iFurlanetto et al.ll2006[ ) in that it char- 
acterizes the transition in terms of the star formation 
rate density, psfr,, a quantity that is also constrained 
by optical/infrared observations of the z '-^ 10 Universe. 
One goal of i)2.1l is to motivate why it is likely for the 
IGM to appear with Tk < Icmb in 21cm absorption 
soon after the first stars formed (the situation this paper 
primarily investigates). A second goal is to understand 
whether minihalos could source the pre-reionization radi- 
ation background, which is the scenario that leads to the 
differential velocity having its largest impact on the 21cm 
signal. In addition, minihalo scenarios have largely (and 
we argue unjustly) been ignored in previous research on 
the extremely-redshifted 21cm background. Next, §2.21 
provides an introduction to the fluctuating 21cm signal 
and also motivates why the velocity difference between 
the baryons and dark matter could lead to larger fluctu- 
ations than in previous models that had neglected this 
effect. Finally, H2.3\ discusses the detectability of 21cm 
signals. 

2.1. the temporal evolution of this signal 

The 21cm line of atomic hydrogen offers a brightness 
temperature contrast with respect to the CMB bright- 
ness temperature of 

T-- = 41.5.H (1 + 6,) (1 mK,(l) 

/ 1 _ ^CMB \ 

1 Z'^/^ mK, (2) 



= -216a;H(l + <56) 
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where Tcmb{z) is the CMB temperature at redshift z, 
Z20 = {l + z)/2 1, Ts is the spin temperature of the 21cm 
hne (|Fieldlll958D . xn is the neutral hydrogen fraction, and 
6b is the overdensity in baryons. Note that for ensuing ex- 
pressions we will sometimes drop the superscript "21cm" 
in T^^^™. Equation ([2|) was evaluated for Ts equal to the 
kinetic temperature of the gas prior to reheating and af- 
ter thermal decoupling or T^'^ = 0.021 (1 + z)^ K. This 
equation demonstrates that the absolute brightness tem- 
perature of the 21cm signal is likely to peak at times 
when the gas was kinetically cold. 

The spin temperature of the 21cm transition interpo- 
lates betwee n the gas temperature, Tk, and Tcmb{z) as 
(|Fieldlll958n 



(3) 



where Xc describes how well particle collisions couple Ts 
to Tk- This coefhcient is effectively zero in the low- 
density IGM at z < 40 (the redshifts we focus on), and 



1.8 X 10"(l + z)"Va 



(4) 



parametrizes the efficacy of Lya scattering at pumping 
the 21cm transition. This scattering couples Ts to Tk via 
the Wouthuysen-Field mechanism (jPieldl [l958l ) . Equa- 
tion dH) assumes c.g.s. units for J^, where Jq is the av- 
erage photon specific intensity at the frequency o f Lya 
and is given by fe.g.. lPritchard fc Furlanettoll2006[ ) 

Ja = E/-(n) / d/L^_^_^eK,.'), (5) 

5 c{l + z)^ 



47r H{z) 



27 AwHiz) 



(6) 



Here, H{z) is the Hubble expansion rate, e is the 
spatially-averaged comoving specific emissivity, /rcc('T-) 
is the probability that an absorption into the n"' 
Rydberg level of atomic hydrogen - and the re- 
sulting radiative cascad e - produces a Lya photon 
(jPritchard fc Furlanettol I2006[) . Zmax is the maximum 
redshift at which a photon could have been pro- 
duced that was absorbed into the n^^ level at z, and 
v'n = fn(l + + z). The approximation that yields 

equation ^ includes only the n = 2 contribution to the 
summation: Photons that redshift into the Ly/3 reso- 
nance contribute negligibly to the pumping of Ts (these 
photons are destroyed after a few scatterings) and higher 
Lyman resonances than Lya are pumped by photons in a 
more restricted spectral range. Equation ([6]) also makes 
the further simplification that the specific emissivity is 
frequency- independent between 1 and 3 Ry with value e, . 

The evolution of the 21cm brightness temperature 
is thought to be determined by the production of 
ultraviolet photons that pump the 21cm line (affecting 
Xa), X-ray heating of the intergalactic gas (Tk), and 
the reionization of hydrogen (xh). In what follows, we 
discuss the order in which these different processes likely 
impacted the 21cm line: 

Lya pumping: Given equations (|4]), ([6]), and 
the number of 0.75 — 1 Ry photons produced per baryon 



incorporated in stars, Na, we can solve for the comoving 
star formation rate density required to satisfy Xq, = 1 
and thereby couple Ts to Tk: 



[PsfrJ 



1.7 X IQ-'^Z. 



-1/2 
20 



104 



Mq yr"^ Mpc" 



(7) 

Stellar population synthesis calculations find Na « 10^ 
for Po p II stars with a sta ndard initial mass function 
(IMF: iLeitherer et al-lfTOOl and jVg 5000 for a t op- 
heavy Pop III stellar population (|Bromm et al.l[200TI ). 

Let us estimate the epoch at which [psfr]q is sur- 
passed, assuming that the fraction f^, of the baryons in 
halos above mass mh are incorporated into stars such 
that psFR = f*PbdJco\\{mh)ldt. Here, /co1i(to/i) is the 
fraction of matter that has collapsed in halos more mas- 
sive than TTih, for which we use the Sheth-Tormen halo 
mass function. For = {0.1%, 1%, 10%}, the character- 
istic star formation rate, [psfr]q, would have been en- 
countered at {6, 15, 19} if star formation traced the mass 
in halos that could cool atomically. These estimates as- 
sumed that such halos have virial temperatures of Tvir > 

10* K, yielding halo masses of niH > 3 x lO'' Z^^^"^ ■ 

Most previous discussions of 21cm radiation have as- 
sumed that halos that cool atomically are the main driver 
of th e Lya pumping (jGnedin fc Shaveiil2004l : iFurlanettol 
l2006f l. However, stars in smaller "minihalos" (halos 
that can only cool by molecular hydrogen transitions) 
may contribute to the production of ultraviolet pho- 
tons. Section [3] will show that ^ 10^ Mq halos must 
generate much of the 0.75 — 1 Ry photons for the 
differential dark matter-baryon differential velocity to 
modulate the local star formation rate in a detectable 
manner. Minihalos with circular velocities of at least 
3.7 km s~^ are massive enough for the gas to be able 
to cool by mole cular hydrogen (iTtegm ark et al.l 119971 : 
lAbel et all 120021: iMachacek etall [20011 ) . and they cor- 

respond to halo masses of > 3.5 x 10 Z2Q Mq. Such 
halos would have achieved = 1 at z = {13, 22, 28} 
for = {0.1%, 1%, 10%}. If instead the bulk of the star 

3/2 

formation occurred in more robust 1.5 X 10*^^20 ' ^^o, 
6 km s^^ halos, Xa — ^ would have occurred at somewhat 
lower redshifts of z = {11, 20, 25}. 

The Lyman- Werner band background (e.g., the spec- 
tral band 11.2 — 13.6 eV) acts to destroy the molecular 
hydrogen and sterilize star form ation, esp ecially in the 
lowest mass minihalos (Haiman et al.lll99 7. 2000). The 
numerical calculations of Machacck et al.l (|2001l) . which 
were confirmed in Wise & Abel (2003), found that 4% 
of the gas is able to cool (and collapse to much higher 
densities) in halos with mass 



merit = 2.5 X 10^ 



8.7x10^^0^^21 



Mr: 



(8) 



where Flw,21 is the Lyman- Werner intensity integrated 
over solid angle in units of 10^^^ erg s^^ Hz^^ sr~^.^ 
The critical star formation rate to raise merit turns 
out to be much less than [psfrIq: Xa — 1 requires 
a radiation background in this band with intensity of 
^"0,21 = 47r/it'Lya (1 + z)/{1.8 X 10") X 10^1 in the same 

^ O'Shoa & Norman (2008) found a similar relation, but with a 
~ 3 times higher normalization. 
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units as Flw,21 (e.g., eqn. S]). Plugging in the numbers, 
= 1 yields a Lyman- Werner intensity of 

Flw.2i ~ 2OZ20 exp[-rLVK], (9) 

noting that Flw,21 ^ -Fa, 21 and tlw is the typical in- 
tergalactic opacity for Lyman- Werner phot ons which can 
can b e 1 — 2 in the absence of dissociations (jRicotti et al.l 
120011) a nd can be la rger once the first HII regions have 
formed ([Johnson et al. 2007). Once Xa — I and as- 
suming equation (|S]), the Lyman- Werner background 
is sufficient to suppress cooling in halos with masses 
< 4 X 10'^ expf--0.47Tf^ tv1 Mq.^ 

I Wise fc Abell ()2008f ) studied the star formation rate 
density in cosmological simulations that followed Pop III 
star formation in minihalos, including self-consistently 
Lyman- Werner radiation backgrounds. They found 
that star formation rate densities sufficient for Xa — 1 
occurred by 2; = 22 — 25 in their simulations. 

X-ray heating: Another critical juncture in the 
evolution of the 21cm signal occurred when psFR was 
sufficient for X-rays to have heated the gas above 
the CMB temperature. Penetrating X-rays are likely 
the most efficient mechanism for reheating the IGM 
(jChen fc Miralda-Escudel2004l ) . However, relating X-ray 
production to the star formation rate (SFR) is more un- 
certain than relating ultraviolet emission to SFR since X- 
ray production depends on the abundances of X-ray bina- 
ries and supernovae . To do so, we follow the methodology 
taken in iFurlanetTol ()2006( ) , using relations calibrated on 
low-redshift galaxies between X-ray luminosity and the 
SFR. In particular, the critical SFR density to heat the 
IGM with X-rays by an amount of Tcmb{z) is 

,s™Jx^4.0xlO-Z,^/^(^)-^(g)-^ (10) 

X ( Lx/SFR_\-' 

V 1040 erg s-i yry °^ 

where fx is the fraction of e nergy that heats the IGM 
(|ShuU fc van Steenberd [19851) . Lx/SFR is the ~ 0.1 - 
2 keV luminosity per unit SFR, tsFR is the timescale 
over which the emitting population had been active, and 
tH = H{z)~^. Equation ^ evaluated Lx/SFR at 
10^0 erg s~^ Mq~^ yr, which is a factor of ~ 5 higher 
than low redshift measurements fo r the same relation- 
ship except between 2 — 10 keV ([Grimm et al.l 120031 : 
IMineo et al.ll201l[ ).^ Interestingly, there is no evidence 
for ev olution in the Lx/SFR, even to z w 6 ([Gowie et al.l 
1201 Ih . Using the Sheth-Tormen mass function, the e- 
folding time, [dlog fco\\{mh) / dt]'^^ , is 0.11 and 0.06 of 

® In addition, if an appreciable number of solar mass stars form, 
these produce infrared radiation that dissocia tes H~, preventing 
the formation of molecular hydrogen. iWolcott-Green & HaimanI 
{2012) found that such an infrared background becomes more ef- 
fective at preventing star formation than the Lyman- Werner back- 
ground if ~solar mass stars comprise > 90% of the stellar mass. 
This would require the high-redshift IMF to be more bottom heavy 
than the 2: = IMF. 

The spectral index of the X-ray emission is uncertain, but 
empirical determinations at low-redshifts are con sis tent with hav- 
ing equal e nergy per log in frequency CRephaeli et al.l 119951 : 
ISwartz et al. 2004). Low-redshift X-ray emission that traces star 
formation is dominated by high-mass X-ray binaries. 



a Hubble time at z = 20 for — 10® Mq and 
rrih = 10^ M0, respectively. 

Is it possible for minihalos to have also dominated the 
X-ray reheating of the Universe? For our fiducial param- 
eters, Flw,21 would have been even larger by a factor 
of 20 when [psfrIx was satisfied than when [psfrIq 
was. Larger fx will reduce the resulting Fi,w,2i, making 
it more difficult for the Lyman- Werner background to 
sterilize minihalos. However, X-rays can also catalyze 
the formation of H2, combatting its destruction via the 
Lyman- Werner background (Haiman et al. 1997). We 
find that at the epoch when [/9sfr]x is satisfied the 
electron fraction is increased by a factor of 10 over the 
relic electron fraction at the cosmic mean density. A 
10 times higher electron fraction means that 10 times 
more molecular hydrogen will form, such that a 10 
times larger Lyman- Werner background is required to 
yield the same H2 fraction. However, this estimate 
represents an upper bound on how much the critical 
Flw,21 of a halo would be increased by X-rays, as the 
recombination time is ^ 1 Hubble time for 5h = 200 
gas with the relic electron fraction at z — 20. Thus, it 
is difficult to significantly boost the electron fraction 
in dense, star f o rming regions with X-rays. In fact, 
iMachacek et al.l ([20031 ) found in simulations that an 
X-ray background only mildly promotes the formation 
of molecular hydrogen in halos. 

Ionizations: The final effect that stars have on the 
21cm signal is via their ionizations. If a stellar popula- 
tion produces A^ion hydrogen ionizing photons per stellar 
baryon, the critical SFR density required to reionize the 
Universe to an ionized fraction of xi is 

[psFR]io„ = 4.4 X 10-1 X, (^) ' 

X ( l^^] ' Mq yr-i Mpc-3, (11) 
V 0.1 4000y °^ ^ ^ ' 

where f^sc is the fraction of ionizing photons that escape 
from their sites of productio n into the IGM. The factor 
/esc i s highly uncertain (^e.g.. iKuhlen fc Faucher-Giguerd 
I2OI2D and likely to be <C 1. For Pop II stars with 
a Scalo initial mass function (IMF), A^inn ~ 4000 for 
Z = 0.05 Zq (Barkan a fc Loebl 120011 ). This number 
varies at the factor of 2-level when changing assump- 
tions regarding the metallicity and the IMF, at least for 
empirically-determined IMFs. However, Pop HI stars 
with a top-heavy IMF are much more effic ient producers 
of ion izing photons, with A^ion ~ 40,000 ([Bromm et al.l 
[2OOI . 

Thus, we find [psfrIq < [psfrIx < fPsFR.lion i n agree- 
ment with the progression found in iFurlanettq (^006) . 
With this ordering, radiation from star formation first 
coupled the spin temperature to the gas temperature 
such that the 21cm signal would appear in absorption. 
Next, radiation associated with star formation reheated 
the Universe, and, lastly, ultraviolet r adiation from stars 
reionized the cosmic gas. In fact, iPritchard fc Loebl 
([2OIOD demonstrated that one could change both Na and 
fx by orders of magnitude and still have a phase in which 
the 21cm signal was in absorption owing to Lya coupling 
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Fig. 1. — Model history of the gas temperature and spin temper- 
ature (top panel) and of the mean 21cm brightness temperature 
(bottom panel) for the parameters Na = 10**, fx = li /* = 0.02, 
Mon = 4000, and /esc = 0.1, assuming that star formation traces 
the mass in atomic cooling halos. The shaded regions qualitatively 
delineate the phases where different radiation backgrounds drive 
the signal: first radiative pumping by far ultraviolet photons, then 
heating by penetrating X-ray ionizations, and lastly reionization by 
extreme ultraviolet photons. The LED A and DARE instruments 
aim to constrain this signal between 10 < 2 < 30. 

Tg to Tk at times prior to X-ray reheating above Tcmb- 
The amplitude of the 21cm signal is larger during this 
phase compared to its amplitude at other times (eqn. [2]). 
The discussion in this paper focuses on this phase. 

Figure [T] shows the history evolution of Tk, Ts and the 
21cm brightness temperature for the parameters /x = 1, 
/* = 0.02, and /esc = 0.1 (defined in eqn.s[2l[10l and[TTl), 
and assuming that star formation traces atomic cool- 
ing halos. The shaded regions qualitatively delineate 
the phases where different processes drive the signal. In 
this model, the 21cm signal appears in absorption until 
z « 11.5 (bottom panel, Fig. [IJ. The temperature his- 
tory for a model in which minihalos dominate the back- 
grounds can be similar, but the minihalo case is more 
likely to be important as redshift increases. 

Our discussion has ignored shock heating as a con- 
tribution to the gas temperature. The Universe could 
heat owing to structure formation shocks. Such heating 
would also suppress the absorption signal that occurs at 
z « 12 — 22 in Figure [TJ Shock heating would be un- 
fortunate as this dip (the sharpest and strongest feature 
that is expected in the sky-averaged 21cm signal) is being 
targeted by the 21cm instruments LEDA and DARE. Its 
presence is also important for the effects studied in this 
paper. In the only previous numer ical study of this dip's 
existence, IGnedin fc Shaveil (|20M ) found that it was sig- 
nificantly impacted by shock heating. Fortunately, our 
simulations ( which are better suited f or this calculation 
than those in IGnedin fc Shaveil I2004D show that shock 
heating does not significantly impact the temperature of 
the intergalactic gas at these redshifts (see Appendix E| . 

Relating the critical psFR derived in this section to 



the psFR measured in rest-frame ultraviolet observations 
is also helpful for gauging when these difi'erent Psfr 
thresholds are satisfied. In particular, iBouwens et al.l 
(|2010[) found psFR « 10"^ Mq Mpc"^ yr"! at z = 8 
as well as that psfr was decreasing quickly with in- 
creasing redshift. A star formation rate density of 
10~^ Mq Mpc~^ yr^^ is only just sufficient to reion- 
ize the Universe in one Hubble time for fcsc = 0.5 - a 
value that is higher than anticipated ~, but it is more 
than sufficient to satisfy [psfr]q- At face value, the 
rest-frame ultraviolet determinations of psFR also sug- 
gest that the psFR thresholds outlined in this section 
are satisfied at relatively low redshifts. However, it 
is l ikely that faint , dwa rf galaxies that are missed by 
the IBouwens et ahl ()2010D observations contribute signif- 
icant l y to the true value of psfr. a.t z > 8 (Bouwens et 
l20lltlKuhlen fc Faucher-Giguerell20ia) . 

2.2. Fluctuations 

Interferometric efforts targeting redshifted 21cm radi- 
ation are not directly sensitive to the mean 21cm signal, 
but instead to spatial fiuctuations in the 21cm brightness 
temperature. In the limit of small fluctuations (appropri- 
ate for the models we consider), the observed brightness 
temperature contra st with respect to Tcm b can be ap- 
proximated as (e.g.. iFurlanetto et al.ll2006l ) 

^j.21cm ^ y,21cm ( I + Sb + -^S^ 
V l+Xa 

+ ^"^"7 Sr-S.X (12) 

where T^^icm jg ^j^g spatial average of equation ([T]) and 
bars over other quantities also denote a spatial average. 
Equation (fT^ is valid at z < 40, once collisions can no 
longer pump the 21cm line, and prior to reionization. 
Respectively, Sb, Sa, 6t, and S\/y are the overdensities 
in baryons, Lya radiation, temperature, and the line-of- 
sight derivative of the line-of-sight proper velocity (aris- 
ing from redshift-space distortions). 

At lowest order, the overdensities in the temperature 
and in the intensity at the Lya resonance can be re- 
lated to the overdensity in baryons and the square of the 
baryon-dark matter velocity difference, w^^, as 

da^{bs,aSb + by2^Sy2)-kWa, (13) 

ST^{bs,TSb + by2^TSv'^) , (14) 

where 

6,2 ^ [vUx)/al - 1] . (15) 

I Dalai et al.l (|201Cl[ ) demonstrated that the expansion to 
first order in 6,2 is a good approximation if 6a and 6t 
are local functions of It'bcl, at least on scales where fluc- 
tuations in v^^ are much less than unity. In addition, 
bs,a is the bias of star forming regions, bs^T is the bias of 
density-tracing temperature fluctuations (equal to 2/3 
for adiabatic evolution)® and by2 a [bv2 ^] is the bias 
with which star formation [heating] traces 6^2 . Lastly, 

* Both temperature "biases" formally enter in convolution at 
times when X-ray heating was important, reflecting the propaga- 
tion of X-ray photons. 
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cr^j. = (Wbc}«: where (...)„ denotes an ensemble aver- 
age over the Maxweh-Boltzmann probabihty distribution 
of Vhc- In addition, Wa{r) describes how Lya photons 
travel from a point source to a distance r away and is 
normalized to have unit norm, entering equation (jl3p in 
convolution. It is given by a simila. r expression to equa- 
tion © (see lBarkana fc Loeb|[2005l) . 

Interferometric 21cm fluctuation efforts aim to mea- 
sure the power spectrum of the 21cm brightness temper- 
ature fluctuations, P2i{k) — (|(5Tf,(A;)p), where tildes sig- 
nify the Fourier dual.^ Combining equations p^ . ([T^ . 
and (HI]) yields 



1 + x„ 



T, 



CMB 0(5,T 



Tk-T, 



CMB 



Wa{k) TcmbK^,t 



1 



Tk-T, 



CMB 



Ps{k) 



P„2(fc), (16) 



where Pi = {\6{k)\^) andP„2 = (|5„2(/c)|2). We are work- 
ing in the limit in which the baryons trace the matter 
such that 6i, — S, and we have dropped the much smaller 
terms that are cx {S{k)'^Sy2(k)), even though they for- 
mally are the same order as P,2 The terms with 
owe to redshift-space distortions UKaiseri [19871) . where 
fi = n ■ k/k and n is a unit vector that points along the 
line-of-sight. 

Previous analyses aside from iDalal et all (|2010l ) had 
not included the P^2 contribution to P21. This new 
term's amplitude peaks at smaller k than P5. In fact, 
even if the coefficient that couples it to P21 is 10'* times 
smaller than the analogous coefficient for Pg, it still 
would contribute comparable power to this other term 
at fc = 0.1 Mpc~* - roughly th e scale at which 21cm 
experiments are m ost sensitive (|McQuinn et al.l 120061 : 
iParsons et al.l 120111) . The Fourier transform of P„2 is 
equal to 



{6y2{x)6y2{X + r)) = gV-lW^ 



2 



where 



ipiir)-: 



2 



3 
3 

^bc "'0 



k^dk 
k^dk 



Pv{k) 



jijkr) 
kr 



Pv{k)j2{kr), 



Pv 



(l^bc(fc)P). 
Wbc(fe) - 
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and 
. a k 



27r2 



Tc{k, a) - Tb{k, a) 



(17) 



(18) 



(19) 



(20) 



(21) 



A dot over a variable signifies differentiation with re- 
spect to time. Tc and Tf, are the baryonic and dark 

^ Any significant contribution to the power spectrum from d^c 
also means the signal is highly non-Gaussian such that higher- 
order statistics such as the trispectrum w ill be easier to detect than 
previously anticipated IIYoo et al.ll201ll) . However, we concentrate 
on the power spectrum here. 



matter transfer functions that map from the primordial 
overdensity^ri, to the overdensity in these components. 
Equation (|17|) can be derived from equation (|2T|) using 
Wick's theorem and the correlation function between the 
cartesian components of the velocity field: 



2 

{v,{x)v,ix + r))„ = ^ [Mr) Sf, + Mr) 



(22) 

where is the Kronecker delta. 

One of the primary objectives of this paper is to calcu- 
late the contribution of the term proportional to Pu2 to 
P21, namely to calculate by2 „ and &„2 ^ in equation (|1( 



An estimate for &i,2 „: If star formation traces the 
amount of matte r that collapses int o halos above some 
mass cutoff, fr^ii. lDalal et al.l (|2010( ) found that 



(^'bc/coll)'t' 
(/cou)d 



1 



'he 



(23) 



provided a good fit to the large-scale bias found in a full 
calculation for the source clustering in the model they 
considered. An alternative bias coefficient comes from 
equating the variance in /coii as a function of Wbc to ((5^2 
(which equals 2/3), and is given by 



(/coll)t' 
(/coll)? 



- 1 



1/2 



(24) 



This alternative bias should hold to the extent that the 
Wbc^drivcn fiuctuations are linear on all scales. Such a 
condition does not hold for the density field at z < 20, 
but is more valid for 5^2 . 

The exercise of plugging fcowivhc) {fco\i)v (1 + 
61,2^ (5^,2) into the right- hand-side of equations ([53]) and 
(|24|) reveals that both expressions for the bias are exact 
in the limit that higher than linear order terms are 
subdominant. We find that either choice of bias agrees 
within ~ 10% among the models we co n sider. However, 
because the calculations in IDalal et al.l ()2010l ) find that 
the bias given by equation ((23)) fares excellently even 
in the somewhat nonlinear cases they considered, our 
subsequent calculations will use equation (j23p for hy2 ,^ 
and for 6^2 when radiation backgrounds source the 
fiuctuations. 

An estimate for 6„2 We derive 61,2 with a similar 
logic to how equation ([M)) was derived: If temperature- 
overdensity fluctuations were coupled to with root 
mean square (RMS) fluctuation amplitude of ctt when 
averaging over regions, then 



(25) 



in order for the variance of the field by2 rp 8^2 to equal ctt- 
As an example of heating, if the fraction /3 of v^^ were 
locally thermalized in shocks (such that the region was 
heated by the factor 1 -I- 5/37^^/9), 



3^2, T 



5/3 



(Kc)„/ 



1 



5/? 



(-^bc), 



(26) 
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2.3. detectability 

It is not necessarily more difficult to detect the fluc- 
tuating 21cni signal from the era of focus here, z ^ 20, 
than from z w 8 - the epoch most efforts are target- 
ing: The absolute brightness temperature of the signal is 
likely to be much larger at z '--^ 20 than at z ~ 8, which 
compensates for t he increased brightness temperature of 
the foregrounds (jPritchard fc Loebl[2008l ). 

For subsequent estimates of the 21cm fluctuating sig- 
nal, we quantify its potential to be detected by con- 
sidering an interferometric conflguration similar to the 
Murchison Widefield array (MWA) but that is optimized 
to the higher redshifts in question (by dilating the sizes 
of the dipoles and baselines). We also consider an in- 
terferometer with 10 times the collecting area that is a 
scaled up replica of the MWA. See Appendix IB] for addi- 
tional details regarding these hypothetical instruments. 
The statistical sensitivities of these instruments (exclud- 
ing sample variance) to P21 measured in spherical shells 
of width Afc = 0.2 k are shown in the ensuing 21cm 
power spectrum flgures. These calculations assume an 
observing time of 1000 hr and a bandwidth of 5 MHz. 
They also assume that the foregrounds can be removed 
from all wavevectors that do not utilize the zero mode 
along the line-of-sight such that the sensitivity at these 
wave vectors is limited by thermal noise (McQuin n et al.l 
|2006( ). For B = 5 MHz, an instrument can at best ob- 
serve k > 0.05 Mpc~^, corresponding to the fundamental 
line-of-sight wavenumber for the chosen bandwidth. 

Larger bandwidths (which allow a measurement at 
smaller k) are desirable to probe the k where Py2 has its 
largest impact on the P2i{k). However, the bandwidth 
(redshift) interval must also be chosen so that there is 
not signiflcant temporal evolution in the signal across 
the band. The amount of collapsed mass evolves on a 
timescale of Az/(1 + z) ^ 0.1 (Section 12. 1|) . which is 
comparable at z = 20 to the Az covered by a bandwidth 
of 5 MHz. Thus, we expect the amount of star forma- 
tion evolves substantially across the redshifts covered in 
this instrumental band. This evolution would result in 
spurious large-scale power at the [T^^'^'^J^-level, which 
dwarfs the signal we are discussing. Thus, extracting the 
Ubc~induced signal from the lowest fc to which a survey is 
sensitive will be challenging (which happen to be where 
the effects of Vhc are most prominent). 

3. THE IMPACT OF Vbc ON THE 21CM SIGNAL 

This section provides analytic estimates for the im- 
print of Wbc on the 21cm signal. These estimates are 
calibrated with the cosmological simulations discussed 
in Paper I. They also require a calculation for how Ubc 
alters the linear growth of density fluctuations. To do 
so, we solve the linear s ystem of equations given in Pa - 
per I, and originally in iTseliakhovich fc Hiratal ()2010D . 
for the growth of modes in the presence of a nonzero 
Vbc- This calculation is initialized with z = 1000 transfer 
functions from the CAMB code.^*' It assumes an instan- 
taneous kinetic decoupling of the baryons f rom the CMB 
at this redshift (as mo tivated in § 5.3.1 in IHul[T995. and 
lEisenstein fc Hui[l998() to solve for the growth of struc- 
ture in two gravitationally coupled fluids, the dark mat- 
ter and the baryons. This approximation for the growth 

^° |http : //camb . Info/ 1 




k [Mpc"'] 

Fig. 2. — Impact of ii^c on the total matter and baryonic power 
spectrum at 2 = 20. Shown are modes in which A4bc ■ ^ equals 
the quoted number. The top panel shows the total matter power 
spectrum, and the bottom panel shows the power spectrum of just 
the baryons. For reference, 95% of space has A4bc > 0-8, 70% 
has A^bc > 1-6, and 5% has A^bc > 3.9. Also shown in the top 
panel is the squared Fourier transform of the top-hat in real-space 
window function (normalized to 0.1 at = 0) that in the spherical 
collapse model would collapse into halos with circular velocities of 
10 km s~^ and 3.7 km s~ - the thresholds at which atomic and 
molecular hydrogen cooling become important, respectively. 

of modes avoids the full Boltzmann code calculation. We 
find that solutions with this approximation excellently 
reproduce the growth of modes in CAMB for the case 
CAMB solves, Wbc = 0. 

Figure [2] features the results of this calculation for 
z = 20, showing the small-scale power spectrum of the 
total matter (top panel) and the baryons (bottom panel) . 
A nonzero Vhc acts as an effective pressure in the baryons, 
decreasing the fc above which pressure smooths fluctua- 
tions and in an anisotropic fashion. The linear variance 
in the baryonic density contrast is altered at the ~ 20% 
level at z = 20 by Vbc- Nonlinear evolution further am- 
plifies these differences. 

3.1. temperature fluctuations from structure formation 
and shocking 

It is possible that shocking or structure formation 
changed the temperature of the Universe between regions 
with different Vhc sufficiently to result in an enhancement 
of the 21cm signal. The period when such heating is most 
likely to source a significant component in the 21cm sig- 
nal is after Lya pumping started to couple the spin tem- 
perature to the temperature of the gas but before the gas 
was reheated above the temperature of the CMB by X- 
rays. Figure [3] illustrates how the 21cm power spectrum, 
ATfc(fc)2 = fc3p2i/(27r2), would be affected if 0.33, 1, and 
3.3% of the energy in Vhc were thermalized using equa- 
tion ()26p and assuming ^ 1. (The relative amplitude 
of the curves do not depend on Xa.) If just a small frac- 
tion of the kinetic energy in Vhc were thermalized, this 
heating would lead to large temperature fluctuations and 
a larger 21cm signal than has been anticipated in models 
that ignore Vhc- We showed in Paper I that the gaseous 
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x„» l,z = 20 



5% Tjff fluctuations / 3.3%of[Vj,J thermalized, P„2 
1.5% Tjf, fluctuations / 1% of [V|,J- therraalized, P,,2 
0.5% Tj|-|- fluctuations / 0.33%) of [v^J" tliermalized, P,2 
MWA-like x 1,10 sensitivity 



10"" 



10- 



k [Mpc J 



Fig. 3. — P„2 and Pg contributions to the azimuthally-averaged 
21cm power spectrum, where ATi,(k)^ = {2tt'^)~^ P21, in cases 
where the effective gas temperature, T^s = (1/Tk)^/, is modu- 
lated by the local value of ?)bc- The different curves represent the 
case where the standard deviation in temperature, crj', is the spec- 
ified value or where the specified fraction of v^^ is thermalized. 
While the curves are for z = 20 and Xa 2> 1, the relative ampli- 
tudes of the P^2 curves are redshift independent, and their overall 
amplitude scales as f^/{fj^ — Tcmb)'^- The amplitude of the Pg 
curve scales with the additional factor of (1 + z)~^. The diagonal 
lines are the statistical sensitivity to the power spectrum in bins 
of Ak = 0.2 A; for an array that has 1 and 10 times the number of 
antenna elements as MWA, assuming an observing time of 1000 hr 
and B = 5 MHz. 



dynamical friction timescale is shorter than the Hubble 
time for dark matter overdensities with masses greater 
than 10'* Mq - roughly halos massive enough to signifi- 
cantly perturb the gas density. The short dynamical fric- 
tion time results in the streaming gas decelerating into 
the dark matter potential wells and in visually apparent 
shocks throughout the simulation volume (see Fig. 4, in 
Paper I). In fact, by z = 20 in the simulation, all gas with 
Sb > 2 has lost its relative velocity owing to dynamical 
friction. The amount of energy dissipated in this decel- 
eration is enough to significantly heat the cosmic gas (by 
a factor of 3, on average) if it were all thermalized. How- 
ever, it is difficult to estimate from first principles the 
amount of energy that this process would actually ther- 
malize: This heating may primarily impact overdense 
gas, which occupies a small fraction of the cosmic vol- 
ume, and the low Mach number shocks sourced by Vhc 



are inefficient entropy generators.** 

Finite Vhc suppresses the collapse of gas onto halos, 
which also impacts the gas temperature. If the suppres- 
sion of structure leads to temperature fiuctuations be- 
tween regions with different Vhc that have standard de- 
viations of ctt = 0.5%, 1.5%, or 5%, each would lead to 
a contribution to P21 given by one of the P^2 curves in 
Figure [3] We can estimate the level of temperature fluc- 
tuations that result from the impact of Vhc on structure 
formation. In the limit that T^^icm ^ p/Tx (that applies 
when Xq, S> 1) and pure adiabatic evolution such that 
Tk oc (1 + 5f,)^/'^, the average 21cm brightness tempera- 
ture depends on Vbc as 



T-(^;bc)=^^f^ (1 + <5,) 



\l/3 



^r2*( 1 



fixed Dbc 



(29) 



where as(vhc) is the standard deviation in the gas den- 
sity in a region with differential velocity Vbc- Thus, 
the average temperature is slightly lower in regions with 
larger variance, at least when as^ <^ 1 so that higher 
order moments are not important. Computing cry — 

{Tb\vbc)/{T^^}v - l)y^ using linear theory yields 0.002 
at z = 20 (a time when the linear theory standard de- 
viation in Sb is ~ 0.5). At z = 20 and for Xa ^ 1, 
the imprint of ctt — 0.002 on the 21cm power spectrum 
would lie a factor of 6 below the 0.5% curve in Figure [3l 
and, thus, be quite small. 

Equation neglects the impact of Compton heating, 
which results in deviations from pure adiabatic evolution. 
Compton heating changes the temperature of a region 
by 0.6 K per Hubble time at z « 20 for the cosmic free 
electron fraction of 2 x 10"'^, although the heating rate 
scales as (1 -I- z)"*. Because the Compton heating rate 
is density independent (noting that the electron fraction 
is largely uniform) whereas the adiabatic cooling rate 
scales as S'^^^ , the fractional temperature in voids is more 
impacted by Compton heating than in filaments. This 
results in a Wbc-dependence to the heating as larger Vhc 
results in shallower void depths. 

Thus, shocks, structure formation and Compton heat- 
ing result in Wbc^dependent temperatures. Cosmological 
simulations are our best recourse for a more accurate 
estimate. We use the cosmological codes GADGETS 

For low Mach number shocks, the energy lost in deceleration 
may not necessarily go into heating the gas. The temperature 
increase across an adiabatic shock is 



To 



16"^ {5M'^ - 1){3M''^ + 1) 



as A1^1 + 



M. 



(27) 



More fundamental is the entropy jump, where we define entropy 
here as S = p/p^i'^. In this case, 



S_ 



-8/3 



{5M^ - 1)(3M-'^ + 1) 



5/3 



as A4— >1 + 



5(M - if 
6 



(28) 

which is zero at first and second order in — 1, where A4 is 
the shock Mach number. For an = 2 shock, {T/To,S/So} 
equals {2.1, 1.2} and, for an = 4, it equals {5.9, 2.6}. Thus, the 
regions where vVtbc is largest may lead to the most heating (even 
though dynamical friction is less effective there). In addition, if the 
simulations fail to capture M — 1 = 0.1 shocks, they are missing 
only 10^'^ jumps in entropy. 
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(iSpringel et alj |2001[ ) and Enzo v2.1.1 (jO'Shea et alj 
l2004fl with initial conditions that include t^bc self- 
consistently (in addition to other improvements; Pa- 
per I). GADGET solves the equations of fluid dynam- 
ics with the smooth particle hydrodynamics method, 
whereas Enzo is a grid code with adaptive mesh refine- 
ment (AMR). Both the GADGET and Enzo codes have 
been shown to conserve entropy at the part in 1000- 
level for the tes t case of the expa nsion of a homoge- 
neous Universe (jO'Shea et alJ 120051 ). This level of en- 
tropy conservation is unusual for hydrodynamics codes, 
and it owes to the entropy-conserving formalism of GAD- 
GET and the 3'''^-order accurate in space, 2"'^~order in 
time Riemann solver employed by Enzo. Thus, both 
codes are well-motivated choices for tracking thermal 
effects in the early Universe, and their much different 
hydrodynamics solvers tests the robustness of the re- 
sults. Each curve in Figure |4] plots the difference 
in (l/rK)j/ between both GADGET and Enzo simu- 
lations between a simulation with A^bc = and one 
with A4hc = 1-9, where (...)m signifies a mass aver- 
age. We refer to this difference as (5(1/Tk)J/. Note that 
f^'^ (x[l + TcMB {Tk~^)m] when Xa Tqub/Tk > 1 such 
that (Tk^^) is the gas temperature-weighting relevant 
to 21cm observations during the period of interest (when 
this signal appears in absorption). 

All of our simulations find that the average tempera- 
ture is larger in the simulations with Aihc = than in 
those with larger A^bc- (Note that on the scale of our 
simulations, < 1 Mpc/h, this differential velocity is a 
uniform wind with a single Mach number, allowing us to 
refer to simulations by their value of At be-) The direc- 
tion of the temperature difference indicates that the sup- 
pression of structure formation and the impact of Comp- 
ton heat- ing are the dominant processes that shape the 
temperature rather than shock heating. The dot-dashed 
curve shows 5{1/Tk)Ii in Enzo between a simulation 
with and one without Compton heating. Compton heat- 
ing contributes more than half of this temperature dif- 
ference. 

With both Enzo and GADGET, the temperatures are 
reasonably converged in resolution and in box size be- 
tween the simulations shown in Figure Hl^'^ We also 
find that Enzo is converged in S{1/Tk)]^ between both 
the adaptive-mesh refinement (AMR) and uni-grid sim- 
ulations. We find that the temperature differences 
are roughly linear in Vhc from comparing with also 
the Albc = 3.8 simulations such that we find ctt ~ 
0.4:S{l/TK}J,f- This allows us to draw a horizontal line in 
Figure 13] that corresponds to the lowest amplitude curve 

In addition, we have put the cosmological codes GADGET 
and Enzo through a battery of tests in order to confirm that our 
results regarding the temperature evolution of the Universe are 
robust. We find that the relative difference in the 21cm intensity- 
weighted temperature in the simulations is robust to: i) the frame 
of reference for the relative velocity of the baryons and dark matter 
on the grid at least for Galilean transforms with boost velocity 
comparable to M^c ii) the maximum time-step size, iii) the grid 
size, iv) the box size, v) the number of particles, vi) the number of 
AMR levels in Enzo, and vii) the chosen hydrodynamics solver in 
Enzo. 

This is not the case if we run GADGET in the standard 
setting rather than with the gravitational softening equal to the gas 
softening as done here, where we find the temperature differences 
are twice as large owing to artificial particle coupling. 
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Fig. 4. — Fractional difference between the 21cm intensity- 
weighted temperature in simulations with A^bc = and in those 
with A4hc = 1.9 (corresponding to iibc = 3 km s~^ at 2 = 100). 
Note that the temperature in the simulations with A^bc = is 
higher, and the fractional temperature difference is similar in Enzo 
to this difference in GADGET. The horizontal line shows the dif- 
ference that leads to the signal given by the lowermost P^2 curve 
in Figure |3] The Enzo simulations here are uni-grid, but we find 
negligible differences when compared with the Enzo AMR calcula- 
tions. 

for the Py2 contribution to P21 in Figure |4l For Wbc- 
sourced temperature fluctuations to contribute a compa- 
rable fraction of the 21cm power as the density sourced 
intensity fluctuations, ax ^ 0.015 must be satisfied (Fig. 
[3]). Thus, the simulations yield temperature differences 
that are a factor of a few too small for the Py2 contribu- 
tion to be an 0{1) contribution to P2i- 

3.2. star formation-sourced fluctuations 

While we found that the coupling of non-astrophysical 
processes to Vhc is insufficient to significantly enhance the 
fluctuations in the 21cm background, the local amount 
of star formation can also couple to fbc and modulate 
associated radiation backgrounds, which in turn affect 
P21- We attempt to model this coupling here. We con- 
sider two disparate models for which halos at z ^ 10 — 20 
dominated the SFR: 

Molecular Cooling Halos: In this case, stars are pri- 
marily formed in minihalos with 10^'^ < Tvir < 
10^ K. These halos cool by exciting molecu- 
lar hydrogen transitions. We use a simplistic 
parametrization to quantify how Vhc could impact 
these halos. In particular, in addition to includ- 
ing the impact of Wbc on the halo mass function, 
we assume that all halos with higher circular ve- 
locities than Vcoo\{z) are able to retain their gas 
and form star s prop ortional to their mass. As in 
iFialkov et all (|2011f ). we use 

l^cooi(z) = {KooLo + [av^zfy^^ (30) 

Previous numerical simulations find T4ooi.o = 
3.7 km s~^ - the approximate threshold to cool via 
H2 when Flw,21 ^ 1 and in the ab sence of bary- 
onic streaming - and a = 4.0 (e.g., IFialkov et al.l 
20111. who used fits to the numerical results of 
Stacv et al.l[20lll and [Greif et al.ll2011bn . We leave 
a as a free parameter that is calibrated with our 
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Fig. 5. — Impact of ti^c on halo dark matter mass function for the 
1 Mpc/h, 2 X 768^ particle GADGET simulations (the largest sim- 
ulations we ran). Also shown is the Sheth-Tormen mass function 
for the same v^,^ calculated for a region that has zero overdensity 
on the box scale, as described in the text, to compare on an equal 
footing. 

simulations, and we use Vcooi,o — 3.7 km s^^ as 
well as Kooi.o = 7.4 km s~^. The latter choice 
roughly corresponds to the minimum circular ve- 
locity that can cool for a Lyman- Werner back- 
ground with Xa — 1 (eqn. [8] evaluated with the 
intensity from eqn. |9]). For reference, 3.7 km s~ 
and 7.4 km s^^ correspond to 4 x 10^ Mq and 
3 X 10^ Mq halos at z 20. 

Atomic Cooling Halos: In this case, halos with Tvir > 
10'^ K dominate the tot a l star formation r ate (e.g., 
IGnedin fc Shaved [200l iFurlanettd l2006l) . It is 
unlikely that Vhc significantly impacts the mass 
threshold that cools and forms stars in this case 
because (1) atomic transitions are a robust coolant 
that depends more weakly on the density of the gas 
and (2) the circular velocity of these halos is much 
larger than Vhc- The impact of Vhc on star forma- 
tion in these halos should derive primarily from its 
effect on these halos' mass function. However, Fig- 
ure [2] illustrates that Vbc has a small impact on the 
matter power spectrum at scales relevant to atomic 
cooling halos, so a large iibc-sourced contribution to 
P21 is not anticipated. 

For both the molecular and atomic cooling cases, our 
calculations assume that star formation is proportional 
to collapsed mass above the minimum mass threshold. 

3.3. numerical estimates for the star formation 
suppression 

Ultimately, our aim in this section is to estimate the 
contribution to the 21cm power spectrum that traces P„2 
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Fig. 6. — Impact of Hbc on mass that can cool within 1 Hubble 
time as a function of halo mass, where each marker shows the gas 
mass that meets our criterion within a virial radius from a halo 
centroid. The blue circles are the {1 Mpc/fc, 2 X 768^ particle} 
simulation with A^bc = Oi a-iid the red stars are the same but for 
M]jc = 1-9. The thin curves with the corresponding color show 
the mean in both simulations. The dashed curve is the cosmic 
baryon fraction times the halo mass, and the short vertical lines 
illustrate the cutoff mass in our simple model for A^bc = and 
A^bc = 1-9 and the specified a. Crudely, the impact of Al^c is 
to shift the minimum halo mass with a ^ A in equation lj30|l . but 
Albc also suppresses star formation in more massive systems and 
in a stochastic manner. 

from star formation-sourced radiation backgrounds. We 
use our cosmological simulations to motivate how the 
star formation rate is suppressed. We first investigate 
Wbc's impact on the halo mass function in numerical sim- 
ulations relative to what one would predict for its im- 
pact on the Sheth-Tormen halo mass function (which 
can be calculated directly from the linear theory mat- 
ter power spectrum). Each panel in Figure [5] shows halo 
mass functions calculated from our two largest GAD- 
GET simulations (1 Mpc/h, 2 x 768'^ particle, one with 
Mhc = and the other with A^bc = 1.9)." These mass 
functions were calculated using the friends-of-friends al- 
gorithm with a linking length of 0.2. For comparison, 
we also show the semi-analytic mass function using the 
Shcth-Tormen theory. We include finite box size effects 
in the Shcth-Tormen mass function as described in the 
ensuing footnote in order to compare on an equal foot- 
ing with the simulation mass functions. Both sets of 
curves are calculated with Albc equal to and 1.9, and 
both simulations use the same random numbers to gen- 
erate the density and velocity fields. While the simu- 
lated mass function is slightly below the Sheth-Tormen 
mass function, particularly at z = 25, the suppression 
of the amplitude of the mass function is comparable to 

In Paper I we showed that particle coupling can impact the 
GADGET simulations when A^bc = and this particle coupling is 
alleviated in the simulations with A^bc > 0- 

To account for finite box size effects 
and [S] we multiply the Sheth-Tormen masi 



n Figures [5] 
function by 

nps{rnh\,J^^^~^^^f^^^^)/nps(rnh\rTrn}, where nps(mh,o-x) is 

the Press-Schechter mass function at mass given (Tx, the RMS 
density contrast in a sph ere of rad ius X. This prescription was 
motivated in iBarkana fc Loebl l|2004l ). 
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Fig. 7. — Same as Figure [6l but for the amount of gas that can 
cool in 0.1 Hubble times. 



that predicted by the semi-analytic mass function cal- 
culation. We will use this semi-analytic mass function 
model in subsequent calculations. 

The other parameter we aim to calibrate from the sim- 
ulations is a. This parameter regulates the minimum 
mass of a halo that is able to cool by molecular transi- 
tions. Figures[5]and[7]show the mass in gas within a virial 
radius of a halo of the specified dark matter mass that 
can cool by molecular hydrogen transitions in 1.0 and 0.1 
Hubble times, respectively, at z — 20. The points show 
the individual halos in the 1 Mpc/h, 2x TGS'^ particle sim- 
ulations, and the curves are the average amount of mass 
in gas that can cool in these simula tions. To calculate the 
cooling time, we use the formula in lTegmark et al.) (|1997[ ) 
under the crude approximation that log(l -|-7Vrec) = 1 to 
calculate the amount of molecular hydrogen, where A'loc 
is the number of recombinations. This equality approx- 
imately holds for gas at the virial density of halos at 
z ~ 20, and it allows us to calculate the cooling rate 
from a single simulation snapshot rather than following 
the density evolution of a fiuid element. 

The amount of mass that can cool trends to zero at 
halo masses less than w 1 — 4 x 10^ Mq (Fig.s [6] and 
[7]). Such a threshold was an ticipated from more de tailed 
first star calculations (e.g., IMachacek et al.ll200lD . and 
larger Aihc shift this turn-over mass to higher values. 
The vertical lines show how much this mass would shift 



have followed the col- 
star 



Most other studies of the impact of Db, 
lapse of g as parcels to much higher densities as a proxy for sta i 
formation I IMaio et al.|[20Tll : IStacv et an[20Tll: IGreif et al.1l2011bl) 

In our opinion, it is not necessarily a disadvantage to use our cool- 
ing criteria on as a proxy for star formation rather than following 
the cooling and condensing gas to much higher densities: Firstly, 
simulations that follow gas to much higher densities do not all agree 
on t he character of star formation in the first halos (e.g., Grcif ct al. 
12011a ;). In addition, feedback processes either from stellar HII re- 
gions within the halo (Alvarez ct al. 2006; Yoshida ct al. 2007) or 
the cosmolo gical Lym an- Werner background (Haiman ct al. 2000; 
[Machacek et al.ll20011) drastically increase the complexity of mod- 
eling star formation at the Cosmic Dawn. However, we have run 
the same simulations in Enzo with molecular hydrogen cooling and 
AMR, and found that our simple estimates did roughly reproduce 
the amount of gas that cooled to much higher densities. 
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Fig. 8. — Ratio of the amount of gas that can cool within the 
specified period and form stars in simulations with A^bc = a^nd 
-^bc = 1-9. Other than the one Enzo curve, all others are calcu- 
lated from simulations run with GADGET. Also shown for compar- 
ison are semi-analytic calculations for this ratio using the molecular 
hydrogen cooling model for the specified a for a box size of 1 Mpc/h 
(thin dot-dashed curve) and 0.5 Mpc//i (thin double dotted curve). 
The simulations are closest to the model for a ^ 4 — 6. 



under our parametrization for the impact on molecular 
hydrogen cooling halos for Vcooi.o = 3.7 km s~^ and the 
specified a. The shift in the critical mass that can form 
stars in the simulations most closely approximates the 
shift seen in the a = 4 case. 

However, the ycooi.o ~ parametrization does not ex- 
plain the halo-to-halo stochasticity in the simulations. 
There are far more halos in the Aihc = 1-9 simulation 
that have zero gas under our cooling criteria than in the 
A^bc = simulation. A small component of this suppres- 
sion results from the halos being slightly less massive in 
the Mhc = 1-9 case (see Fig. [S]), but most originates 
from Ubc's suppression of gas accretion onto these halos. 
In Paper I, we showed that one source of this variation 
from halo to halo results from the orientation of filamen- 
tary fiows onto a halo. 

Figure [8] shows the ratio of the total gas that can cool 
in the simulation boxes with A^bc — 1-9 to their coun- 
terpart with A4hc — 0, under the same cooling criteria 
and using three sets of the GADGET simulations and 
one set of Enzo simulations. Also shown for comparison 
are semi-analytic calculations for this ratio in the cases 
with T4ooi,o — 3.7 km and with a = 2, a = 4 or 
a = 6. The simulation curves are most consistent with 
the a « 4 — 6 curves, but note that the simulations be- 
come less reliable at the highest z owing to the rareness 
of these halos. We have also looked at the suppression of 
■Mhc — 3.8 and reached a similar constraint on a. 

In conclusion, the factor by which star formation is 
suppressed is on average similar for both the 1 and 
0.1 cases in FiguresHEl and[Sl We note that the 0.1 
Hubble times criterion is likely most appropriate given 
the exponential growth with time in the number of col- 
lapsing halos. Also, we find a similar suppression in our 
simulations with Enzo as with GADGET. We conclude 
that a = 4 — 6 is the best match to the suppression in 
the SFR under this simple parametrization. While much 
of our discussion focussed on the Vi^ooi.o = 3.7 km s~^, 
it seems natural to assume a is fixed as different physics 
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changes Kooi.O: such as the presence of a Lyman- Werner 
background. We showed that this was the case if we 
changed the coohng time threshold (compare Fig.s[6]and 
[7|). The ensuing discussion will further motivate this as- 
sumption. 

3.3.1. explanation for SFR suppression 

We have found that the simulations prefer a « 4 — 6 
in equation (1301) , which parameterizes the minimum star 
forming halo. Let us attempt to understand this result 
with two toy models. An estimate for the overdensity of 
gas in a virialized minihalo when Vhc — and ignoring 
cooling is 



where is the temperature of the IGM at virializa- 
tion under the idealization of a purely adiabatic collapse 
and hydrostatic equilibrium (Tegmark et al. 1997).^^ In 
the case of finite Vhc and where its energy thermalizes 
in shocks during the collapse, this expression is altered 
such that TiGM.ad -> TiGM,ad(l + 5MlJ9), resulting in 
6b (X (1 + 5M^^/9)~^/^ . There is some threshold virial 
temperature, Tvir,*, that obtains high enough Si, and Tvir 
to cool within a Hubble time. The molecular-hydrogen 
cooling time is proportional to exp[512 K/T]/nH2- ^ 
crude level, we can ignore the factor exp[512 K/T] be- 
cause T4ooi,o = 3.7 km s~^ corresponds to an 800 K halo 
at z = 20 so that this factor is changing by less than the 
factor of 2 above this velocity scale. With this approxi- 
mation, Tvir,* (and also Vjooi) scales as w (1 -I- 5Ml,^/9) 
in this model, since Sf,. This scaling for Vcooi re- 

sults in a modulation comparable to a = 6 at z = 20 
in the parameterization in equation (j30p and with a de- 
creasing with increasing redshift. (Including the factor 
exp[512 K/T] suppresses the effective a somewhat.) 

An even simpler (but related) model uses that veloc- 
ities scale as (1 -I- Sy^^ in an adiabatically collapsing 
region, at least in the absence of dissipative processes. 
Thus, at halo densities 6 ~ 200, the streaming velocity 
accelerates and becomes ~ 6 Vbc, where Vhc is the velocity 
at the cosmic mean density. The total effective pressure 
at a halo virial radius scales as ^ {V^ooi o + t^ t'bc(2)]^ /3}, 
which has a similar form to the parameterization for V^ooi 
given in equation (pO)) and suggests a « 4. If the second 
term is larger than the first, the ram pressure is larger 
than the thermal pressure (potential depth). One can 
also understand this scaling as the requirement that the 
circular velocity of the halo must be larger than the lo- 
cal streaming velocity 6 Vhc) for the gas to be focussed 
into the potential well, where it can shock and cool. 

3.3.2. Ultraviolet pumping 

Now that we have calibrated our model, finding a « 
4 — 6, we can use it to make predictions for the 21cm 
signal. We first focus on 21cm inhomogeneities resulting 
form ultraviolet radiation from the first stars as this is 
likely to occur earlier than other star formation-sourced 
inhomogeneities (as argued in Section |2.1|) . 

Equation l|3H l is start ing to break down at interesting Tvir: 
When Tvir > 260, equation Jstj yields <5j, > 200 at z = 20 (i.e., the 
densities it returns are higher than the virial density). 
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Fig. 9. — Pj,2 (thick curves) and Ps (thin curves) contributions to 
the azimuthally-averaged 21cm power spectrum for « = 4, where 
ATi,(A;)^ = fc'^P2i/(27r ), for different models of how f^c modulates 
the ultraviolet pumping. The top panel examines the case {z = 20, 
Xa = 1} and the bottom {z = 30, Xa = 0.1}. See the text for 
details regarding the plotted models. The sensitivity curves are 
for an MWA-like instrument but optimized for 2 = 20 and for a 
scaled-up MWA with ten times the collecting area. These curves 
assume the same observing parameters as described in the caption 
of Figure |3] For a = 6, the minihalo curves that are sourced by 
P„2 shift upward by a factor of 2. 

Figure |9] shows the P„2 and Ps contributions to the 
21cm power spectrum, where ATi,{k)'^ = P21 / {2tt^) , 
for different models of how Vhc impacts star formation as- 
suming {z = 20, Xa = 1} in the top panel and {z = 30, 
Xa = 0.1} in the bottom panel. These curves assume 
that the X-ray heating and ionizations from these stars 
has yet to impart significant fluctuations in the 21cm 
background - that the only fluctuation sources are den- 
sity and Xa- Each curve can be rescaled to a different 
Xa by noting that this parameter only changes the am- 
plitude of P21 via the factor T^{1 + Xa)~^- (This is not 
the case at A; > 0.1 Mpc~^ for the Ps contribution, where 
the fluctuations begin to be sourced directly by Sb rather 
than Ja as at smaller wavenumbers.) Our calculations 
also assume that the sources' specific emissivity scales as 
e (X ly" with a = 0, e{z) = exp[— (rfcoiii?)""^ ^^/(l + z)], 
noting that r^r^^j = d log fco\\{mh) / dt, dzj (1 + z) = dt, 
and riinax = 20 (see ij2.1l) .^^ We assume the same 
parametrization for f.{z) when considering X-ray heating 

in gS331 

The thick curves in Figure [S] are the P^i contribution 
to P21 and the thin are this but for Ps- The molec- 
ular cooling halo cases with a = 4 are given by the 

We have ignored that that Lyman- Werner background can 
suppress star formation in regions that produce that are producing 
Lycf. Holzbaucr & ; Furlanettg (i2012i ) found that such fluctuations 
are unlikely to have a large impact on P21 . It would decrease the 
amplitude of both the Pj,2 and P^-coupled 21cm fluctuations aris- 
ing from star formation. 

In addition, we have not included the effects of lensing 
and peculiar velocities on <5a and (5y. These effects enter at 1/b 
IIBarkana &: Loe"EII2005l ) and alter the contribution proportional to 
(5;, at the ~ 10% level (but do not quantitatively effect our results). 
Both effects do not alter the component that traces P„2 . 
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dashed curves (Vcooi.o > 3.7 km s ^) and dotted curves 
(K;ooi,o > 7.4 km s~^). The same calculations but with 
a = 6 shift the P„2 curves in the minihalo-sourced mod- 
els upward by a factor of w 2. The atomic cooling halos 
case is given by the solid curves, and it only appears in 
the top panel since atomic cooling halos would not exist 
in sufficient abundance to pump the 21cm line at z = 30. 
In all cases, the relative impact of Wbc becomes smaller 
with time (decreasing redshift) because (1) Pg is growing 
as (1 + z)^'^, (2) halos at all mass scales are becoming 
less rare such that the impact on the growth of modes is 
less accentuated, and (3) Vhc and, hence, its role as an ef- 
fective pressure for a halo with potential depth ^ Vcooi o 
is becoming smaller with time. 

Thus, we find in both considered minihalo cases that 
a significant component of the power is sourced by Py2 
at k - 10~2 Mpc~^ (Fig. |9l). T his is consistent with 
the findings of iDalal et al.l (|2010( ). although our model 
for the impact of Wbc on star formation is more conser- 
vative. Unfortunately, the power is suppressed at larger 
wavenumbers by the free streaming of ultraviolet pho- 
tons, and 21cm instruments are most able to observe 
k > 0.1 Mpc~^ for reasons detailed in ^2.3\ Thus, we con- 
clude that P^2 does not contribute significantly to P21 on 
observable scales at z = 20. For the z = 30 case in Figure 
[9l the Py2 contribution is beginning to become compa- 
rable to the contribution from at fc ~ 0.1 Mpc~^, 
roughly the largest observable scale. However, observing 
the signal at z = 30 is a more challenging venture. 

We can also quantify how much both (1) VhcS effect on 
the suppression of the growth of dark matter structure 
and (2) its impact on disrupting gas accretion onto mini- 
halos (and modulating the minimum mass that can form 
stars) contribute to the inhomogeneities Vhc sources in 
the 21cm power spectrum. We find that effect (2) tends 
to dominate in our models, sourcing w 80% of the P„2- 
coupled signal at z = 20 in both minihalo models. This 
percentile also holds for our ensuing estimates in tj3.3.3l 

3.3.3. X-ray heating 

We argued in Section 12.11 that Lya pumping with 
^ 1 is likely to have occurred before X-rays signif- 
icantly heated the IGM, but, admittedly, there is much 
uncertainty in the high-redshift production of X-rays. 
However, the fluctuations from X-ray heating are less 
damped than those from the Lya background because of 
the shorter mean free path of the X-ray photons that are 
responsible for such heating. The amount of damping is 
small on all scales we plot for the spectral index of the 
background that is assumed (—1 in intensity per unit fre- 
quency), and we found that this result was generic over 
a range of spectral indices. 

Figure [TU] shows predictions for P21 when the fluctua- 
tions are primarily driven by z 20 X-ray backgrounds. 
For simplicity, we assume ^ 1 so that fluctuations in 
Xa are zero and Tk = (Pad + 7cMB)/2. The amplitude of 
our curves can be rescaled to other Xa and Tk. Values of 
Xa ^ i would imply large Lyman- Werner backgrounds 
that would suppress star formation in the environments 
we consider. However, note that changing a:^ = 2 to 
Xa — 00 has little impact on our predictions. In addi- 
tion, we assumed that the critical velocity that can form 
stars is determined by equation (|30]) with a = 4 for the 
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Fig. 10. — P„2 (thick curves) Ps (thin curves) contributions 
to the azimuthally-averaged 21cm power spectrum, AT[,(fc)^ = 
fc^P2i (fc)/(27r^), when inhomogeneous X-ray heating dominates 
the fluctuations. The Hnes correspond to the same star forma- 
tion model as in Figure |9] with a = A for the minihalo cases. Here 
we assume ^ 1, and that X-rays are responsible for raising 
the temperature to Tk = (Tad -I- 7cmb)/2. The fiducial choices 
maximize the amplitude of the component of the signal from X- 
ray heating. The sensitivity curves for an MWA-like instrument 
are the same as described in Figure [3] 

minihalo curves in this figure. If we instead used a — Q 
(which is also consistent with the modulation observed in 
our simulations), this results in a factor of ss 2 increase 
in the normalization of these curves. 

The fluctuations are larger on observable scales in 
the case in which X-ray heating drives the fluctuations 
compared to the case where the fluctuations originate 
from the ultraviolet background. The P„2 contribution 
in both of the minihalo models in Figure [10] may be 
detectable with an MWA-like instrument. However, 
Figure [10] illustrates that only in the optimistic model 
with Vcooi.o = 3.7 km is the Py2 contribution 
comparable to the Pg contribution and then only at 
fc < 0.1 Mpc~^. Thus, even in this optimistic case in 
which X-rays dominate the fluctuations, the signal is 
not dramatically enhanced. 

Thus, the contribution of Py2 to the 21cm signal is 
largest when X-rays from minihalos reheated the Uni- 
verse. However, in models in which atomic cooling halos 
dominate the z 20 X-ray background, the contribu- 
tion to the 21cm power spectrum that is sourced by Vhc 
is likely to be insignificant. 

4. CONCLUSIONS 

This study presented semi-analytic calculations, which 
were calibrated with large cosmological simulations, 
aimed at understanding the early Universe and its 21cm 
signatures. We focused on z 20, an epoch when fluctu- 
ations in the diffuse 21cm background arc not necessarily 
more difficult to detec t than these fluctu ations from any 
other cosmic era (;e.g.. lPritchard &: Loebl [2008). A detec- 
tion of this signal would provide a window into when the 
first stars formed and into a time before the IGM had 
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been reheated by astrophysical sources. 

We focused on the question of whether the dark 
matter-b aryon sup ersonic differential velocity, Wbc 
(jTscliakh ovich fc H irata 20 l"o[). is likely to enhance th e 
2- 20 21cni signal. In iCTLearv fc McQuinnI (|2012i ). 
we showed that typical values for Vhc have a dramatic 
impact on the morphology of structures on scales of 
10 — 100 comoving kpc. Whether Vhc impacted the 21cm 
signal in an observable manner boils down to whether 
v\yc also affected the — 10 comoving Mpc modes to which 
interferometric 21cm efforts are anticipated to be sen- 
sitive. In fact, we showed that if just 1% fractional 
fluctuations in the gas temperature were correlated with 
the large-scale Vbc flows, this coupling could lead to a sig- 
nificant new component to the signal. We investigated 
three sources of Ubc~coupled 21cm fluctuations relevant 
to the epoch after the first stars turned on but prior to 
reionization: (1) heating from structure formation and its 
associated shocks, (2) ultraviolet pumping of the 21cm 
line from the first stars' emissions, and (3) heating by 
the X-ray background produced by the first supernovae. 
X-ray binaries, and miniquasars.^ 

We found that the first of these sources, Wbc's impact on 
structure formation and shocking, did not contribute a 
significant level of large-scale 21cm fluctuations. In par- 
ticular, we found that simulations both with and without 
Ubc did not yield large enough differences in the aver- 
age gas temperatures to result in fluctuations that were 
comparable to the fluctuations in the 21cm signal from 
density fluctuations. We estimated that shocking con- 
tributes < 10% of the power at all wavenumbers. This 
finding held true despite the fact that most of the over- 
dense gas was decelerated into the potential wells of the 
dark matter via dynamical friction in our simulations, 
losing its relative velocity by z — 20. This process leads 
to supersonic wakes and shocking throughout the cosmic 
volume, but not to significant heating. 

The other mechanism by which Wbc could imprint new 
fluctuations in the 21cm background involves the spatial 
modulation by Ubc of the formation of the first stars and 
their associated radiation backgrounds. We investigated 
with a suite of cosmological simulations whether it is 
plausible that Vbc modulates early star formation. At 
z = 20, we found « 3 times more mass in gas that could 
cool in 1.0 (or 0.1) Hubble times and form stars in our 
simulations with A4hc — than in those with Aihc = 1-9. 
Surprisingly, even the amount of dense gas in some of our 
most massive simulated halos (10^ — 10^ Mq) could be 
suppressed at the order-unity level for some halos in the 
simulations with A^bc = 1-9. We provided simplistic 
analytic estimates that reproduced the average amount 

We did not consider the impact of f^c on the signal from the 
reionization epoch as the consensus is that this epoch was driven by 
star formation in the atomic cooling halos that are too massive to 
be im pacted by v^,^ (but see Dalai ct al. 2010 and Bittncr & Loebj 
120111) . vy^^ would also impact the dumpiness of gas and, thus, the 
number of recombinations during reionization (Joanne Cohn, pri- 
vate conversation). This suppression would delay reionization in 
regions with smaller vy^^. However, any X-ray preheating prior to 
reionization would act to eliminate the dumpiness on scales im- 
pacted by fbc, and it is thought that ~ 10^ Mq minihalos — which 
are only moderately impacted by ti^c - are likely to dominate the 
number of a bsorptions from minihalos in all reionization scenarios 
miiev et al.i r2005: McQ uinn et al.l|2007l ). Nevertheless, an investi- 
gation of this coupling mechanism would be interesting. 



of suppression in star-forming gas that was found in the 
simulations. These estimates revealed why the impact is 
significant on halos with circular velocities of « 10 vy^c, in 
contrast to previous estimates. In detail, the amount of 
suppression in the simulations varied significantly from 
halo to halo at fixed halo mass and depended on, for 
example, the orientation of filamentary accretion flows 
with respect to the baryonic wind (see Paper I). 

The fluctuations in the 21cm background from inho- 
mogeneous ultraviolet pumping are likely to be present 
at higher redshifts than fluctuations from other stellar- 
mediated radiation backgrounds. Unfortunately, the 
long mean free path of ultraviolet photons to redshift into 
the Lya resonance dampens the fluctuations from Vbc on 
observable scales. We found that the Ubc-dependence of 
ultraviolet pumping was a significant contribution to the 
21cm anisotropy at k — 0.01 Mpc~^, alt hough not as 
large as in the models of IDalal et al.l (|2010f ) in which the 
impact of Ubc(-z) on the star formation rate was rather 
extreme. However, we found that Wbc could not source a 
significant component of the anisotropy on the largest ob- 
servable scales (fc > 0.1 Mpc"^) at z = 20. At z = 30 (a 
redshift from which the diffuse 21cm background would 
be more difficult to observe), we found that the Vhc con- 
tribution at fc ~ 0.1 Mpc""'^ could be more significant. 

We concluded that Vhc is most likely to leave an 0(1) 
imprint on the 21cm signal if X-ray heating from star- 
forming minihalos drives fluctuations in the intergalactic 
gas temperature. Such X-ray production at high red- 
shifts is very uncertain, and it is by no means guaranteed 
that X-ray reheating was sourced primarily by minihalos 
rather than by more massive halos that could have cooled 
atomically. In fact, most previo us models for the 21cm 
signal a ssumed the latter (e.g. iGnedin &: Shaveil 120041 : 
lFurlanettQ.2006, ). X -ray reheating by minihalos would 
be more likely if minihalos produced more X-rays per 
unit star formation rate at z = 20 than z = galaxies 
(such that X-ray reheating occurred before the Lyman- 
Werner background quenched star formation in miniha- 
los) . Previous studies have made arguments for why 
this could have been the case (IOh|[200ll: iMirabel et al.l 
I20TI . In the cases we considered in which > 3 x 10^ Mq 
and > 3 X 10^ Mq halos formed stars in proportion to 
their mass, we found that the 21cm anisotropy from Ubc's 
modulation of the X-ray background could be compara- 
ble to the 21cm signal sourced by density fluctuations 
at z ~ 20. Most of this anisotropy was sourced by the 
impact of i^bc on the gas accreted by these halos and 
not from the suppression of the halo m ass function (the 
effect considered in Tseliak hovich fc H irata 2010). The 
acoustic oscillations in this Dbc~sourced anisotropy make 
this finding particularly interesting as it results in a more 
distinctive signature for 21cm observatories to target - 
a signature which would indicate that the Universe was 
reheated by star-forming minihalos. 

Finally, we also investigated whether structure 
formation-initiated shocks reheated the Universe in 
Appendix A, which would reduce the absolute brightness 
temperature of the 21cm signal from z ^ 20, times 
when this line is anticipated to appear in absorption. 
The one previous study that investigated this signal 
with cosmological simulations found that such shocks 
would dramatically suppress the amount of absorption 
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(|Gnedin fc Shaved [200l . The recently-funded LEDA 
and recently-proposed DARE 21cni observatories aim to 
detect this 21cni absorption trough in the sky-averaged 
21cm signal, and the sensitivity of these efforts scales 
with the depth of the absorption trough. Fortu- 
nately, we found that shock heating only suppresses the 
21cm absorption at a fractional level of < 20% at z > 10. 
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Fig. 11. — Evolution of the gas temperature in GADGET (thick curves) and Enzo (thin curves; multiplied by 0.6 in the top panel) for 
simulations with the specified comoving box sizes and gas particle/grid dimensions. The top panel shows (l/Tk)^^ , which is how the 
brightness temperature of the 21cm signal weights the temperature when ^ 1. Also shown in this panel is the gas temperature for a 
homogeneous, expanding Universe normalized by eye (double dotted curves labelled "pure adiabatic evolution"). The bottom panel shows 
the fractional residuals of (1/Tk)^j^ with respect to the adiabatic evolution case. Both GADGET and Enzo have almost identical residuals 
when compared at the same resolution and box size. The curve labelled "no Compton" is from a simulation that did not include Compton 
heating but that was initialized with a temperature so that the simulated gas cooled to 10 K at 2 = 20. 

APPENDIX 

A. WHEN DOES SHOCK HEATING RESULT IN THE INTERGALACTIC TEMPERATURE DEPARTING FROM ADIABATIC 

EVOLUTION? 

Here we quantify the importance of shock heating in determining the temperature history of the gas during the 
cosmic Dark Ages. It is plausible that prior to X-ray reheating the Universe shock heated to significantly higher 
temperatures than if it had simply cooled adiabatically with the expansion of the Universe. This heating would be 
analogous to what cosmological simulat ions find at z 5. 1 1 whe re ~ 50% of the simulated mass shock heated to > 10^ K 
owing to 100 km s~^ convergent flows (|Cen fc OstrikeHll999D . However, at z ~ 20, a 0.3 km s^'^ flow was supersonic, 
so weak flows could have sourced shocks and heated the gas. 

The middle and bottom panels in Figure [TT] quantify how much heating occurs in our simulations. The top panel 
shows (Tk"^)^/, where {...)m signifies a mass average, and the bottom panel shows the deviation of (Tk^^)^/ from adi- 
abatic evolution. Note that T^^ oc [1 +Tcmb {Tk^^)m] for Xa Tcmb/Tk ^ 1. Thus, (Tk^^)^/ is the gas temperature- 
weighting most relevant to 21cm observations (after ultraviolet pumping becomes effective). Both GADGET and Enzo 
have almost identical residuals when compared at the same resolution and box size. 

Both codes predict w 10 — 15% deviations above adiabatic evolution by z = 15. About half of this deviation owes 
to Compton heating rather than shocking. (The "no Compton" curve does not include such heating.) We have also 
run several adaptive simulations with Enzo that result in much higher resolution in overdense gas and find good 
convergence with the temperature evolution in these uni-grid Enzo calculations. We conclude that structure formation 
shocks will not qualitatively change the 21cm signal, even in scenarios where the IGM is not impacted by astrophysics 
at z w 10. 

Two prior studies ha d looked at the importance of shock heating in setting the temperature of the early Universe. 
iFurlanetto et al.l (|2006[ ) estimated the fraction of mass that shock heated above the CMB to be /sh ~ {0.1%, 3%, 25%} 
at z = {30, 20, 10} based on the amount of mass that was at turnaround in potential wells that were massive enough 
(i.e., had large enough characteristic velocities) to heat the IGM above the CMB temperature. However, shocks could 
occur even prior to turnaround, as only 0.1 — 1 km s~^ flows can be supersonic and shock. For example, as shells 
cross in voids, this crossing will generate shocks in the gas (e.g., Bertschingcr 1985; although, shell crossing in spherical 
top-hat voids only occurs at a linea r overdensity of —2.7). In support of the possibility that shocks do not just occur 
at or after turnaround, the study of iGnedin fc Shaveil (|2004D conclud ed using cosmolo g ical si mulations that structure 
formation shocks are more prominent than the simple estima t es in IFurlanetto et al.l 1)20061 ) suggest, creating order 
u nity differences even at the highest redshift iGnedin fc ShaveJ (|2004') c onsidered, z = 17. (Although, the simulations 
in IGne din fc Shaver 2004 did not resolve the Jeans' scale.) The IFurlanetto et al.l (|2006D semi-analytic estimate is more 
consistent with the heating seen in our simulations. 

B. SENSITIVITY TO 21CM SIGNAL 

The principle difficulty with detecting high-redshift 21cm radiation is that the sky is much brighter than the 21cm 
signal owing to foreground emission, the dominant foreground being synchrotron from the Galaxy. Foreground emission 
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not only must be subtracted off to isolate the 21cm signal, but also is likely to be the dominant source of statistical 
noise in a measurement of P2i- The anticipated spectral smoothness should allow all significant extraterrestrial 
foregrounds to be separated from the 21cm signal. However, since the brightness temperature of the sky scales as 
Tsky « 240 K (I//150 MHz)"^-^^ (|R,ogers fc Bowman! [2008) . it is thought that detecting the 21cm si gnal will be an 
even more difficult task as the targeted redshift increases (the observed frequency decreases; e.g., IMcQuinn et al.l 
l2006f l. The challenge with statistically detecting the 21cm power spectrum at wavenumber k scales proportionally to 

To quantify the detectability of the signal in our models, we consider two hypothetical interferometer designs in this 
paper. Our first design is motivated by the design of the MWA instrument, which is targeting the 21cm signal from 
z ^ 10. Specifically, for ^ = 20 we assume an instrument that consists of a 40 meter core of closely packed antennae 
with filling fraction 0.5 and with an baseline distribution to larger radii. We assume each dipole contributes A^/4 
in collecting area and that the dipoles have spacing A/2. In addition, we contract or dilate this z = 20 array in 
proportion to the targeted wavelength when we consider the signal from z = 15 or z = 30 in order to optimize the 
instrument to the specified redshift. 

MWA's cost is ~ 20 million US dollars for 500 antennae. The cost of our analogous lower frequency instrument 
would be comparable. We also consider in the body of this paper a next generation array with 10 times the number 
of antennae elements and, for z = 20, a 160 meter core prior to an r"^ falloff, which we denote as "M WA-like x 10". 
The sensitivity of these two hypothetical instruments to the signal is discussed in Section 12.31 See IMcQuinn et al.l 
(|2^006.) for additional details regarding our sensitivity calculations. 



